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Discrete breathers with purely anharmonic short-range interaction potentials localize super- 
exponenttally becoming compact-like. We analyze their spatial localization properties and their 
dynamical stability. Several branches of solutions are identified. One of them connects to the well- 
known Page and Sievers-Takeno lattice modes, another one connects with the compacton solutions of 
Rosenau. The absence of linear dispersion allows for extremely long-lived time-quasiperiodic local- 
ized excitations. Adding long-range anharmonic interactions leads to an extreme case of competition 
between length scales defining the spatial breather localization. We show that short- and long-range 
interaction terms competition results in the appearance of several characteristic cross-over lengths 
and essentially breaks the concept of compactness of the corresponding discrete breathers. 



PACS numbers: 63.20.Pw, 63.20.Ry, 05.45.-a 



I. INTRODUCTION 



Energy localization due to a nonlinearity in dynami- 
cal systems has been observed for more than one cen- 
tury 1], and the effect of an exact balance between non- 
linearity and linear dispersion of wave packets leading 
to the appearance of soliton excitations has become a 
paradigmatic example in nonlinear science which can be 
found in various textbooks. In the past decade remark- 
able achievements in the study of localized nonlinear ex- 
citations were made with the discovery of stable local- 
ized modes in spatially discrete translationally invariant 
Hamiltonian systems - discrete breathers (DBs) . They 
have been proved to be generic exact time-periodic so- 
lutions of the corresponding coupled nonlinear ordinary 
differential equations, eventhough the latter are gener- 
ally non-integrable. It is worth mentioning that discrete 
breathers have been observed experimentally in various 
physical systems including coupled optical waveguide ar- 
rays H, ^1^, coupled Josephson junctions microme- 
chanical cantilever systems 7, 8, 9], anti- ferromagnetic 
crystals jHllll, high- Tc superconductors Discrete 
breathers are predicted also to exist in the dynamics of 
dusty plasma crystals [Tsj . 

Among the most important characteristics of a local- 
ized excitation are its localization length and the spatial 
decay characteristics of its amplitude. Although DBs can 
be localized practically on a single site, in most of the 
cases they have exponentially decaying tails (similar to 
their continuum counterparts - solitons). This is true 
if the interaction potentials are reasonably short ranged 
(see for details). However, when an anharmonic in- 
teraction between adjacent sites is much stronger then 
the harmonic one, localized excitations can become even 
more compact. As it was demonstrated by Rosenau and 
Hyman [l5llT6l |. in continuous systems nonlinear localized 
excitations may compactify, i.e. gain strictly zero tails, 
under nonlinear dis pers ion. The same was conjectured 
for discrete systems [l^ , however later it was shown that 
in discrete systems localized excitations cannot have an 



exact compact structure JJJ , but the tail decay follows a 
super- exponentional law g-t^expfm^ provided that the in- 
teraction is purely short-range. This fact was then con- 
firmed numerically 19], and the corresponding breather 
solutions were coined compact-like 19] or almost- compact 
20] DBs. 

If compact-like breathers are dynamically stable, we 
may expect that localized perturbations of such solutions 
will lead to a quasiperiodic in time evolution, which will 
not induce a radiation of energy away from the breather. 
This is in contrast to the well-known existence of such a 
radiation for systems with linear dispersion |2H . There 
it appears due to the resonant overlap of combination 
frequencies of the internal perturbed breather dynamics 
with the spectrum of small amplitude plane waves. In 
the case of purely nonlinear dispersion the width of this 
spectrum is zero, and thus the origin of the radiation 
is removed. One expects then that perturbed compact- 
like breathers will not radiate energy away, giving rise to 
genuine quasiperiodic compact-like breathers. 

Another important issue which might drastically 
change the rate of spatial decay in DBs tails is the pres- 
ence of long-range interactions, essential e.g. in sys- 
tems with weakly screened Coulomb interaction such as 
ionic crystals, or various biomolecules. Usually decaying 
slower than exponentially in space, long-range interac- 
tions introduce a crossover length as a result of com- 
petition of the two essentially different length scales 14|. 
They can also lead to the appearance of energy thresholds 
for DBs in some cases, where a pure short-ran ge in terac- 
tion would not be capable of producing any. In it was 
demonstrated, that the effect of length-scale competition 
with long-range algebraically and exponentially decaying 
interaction can lead to a new type of multistability of 
DBs, when in a certain model parameter regime several 
different types of DBs coexist having the same value of 
the spectral parameter (i.e. velocity or frequency). 

It is the purpose of this paper to address the 
abovelisted issues. The paper is organized as follows. In 
Sectioniniwe introduce the model. We demonstrate, that 
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the specifically chosen nonlinear potentials allow one to 
completely separate temporal and spatial dependencies 
and thus significantly simplify the analysis of the prob- 
lem. We derive the nonlinear coupled algebraic equa- 
tions for the spatial profile of a solution and in addition 
an ordinary differential equation (Dufhng equation) for 
the master function describing uniform oscillations of all 
the sites with time. In Section IIIII we obtain the dif- 
ferent types of discrete breather solutions. We demon- 
strate, that in general the model supports two classes of 
discrete breathers with completely different dynamical 
properties. We then study the linear stability properties 
of basic types of DBs and observe quasiperiodic local- 
ized excitations. In Section IIVI we reveal the effect of 
long-range interactions along the chain on properties of 
DB solutions. We show that the presence of non-local 
dispersive terms result in the appearance of several char- 
acteristic cross-over lengths. We derive estimations for 
these cross-over lengths, as well as asymptotes for am- 
plitude distribution in DB tails, on the basis of a simple 
three-site model. In Section |3 we conclude. 



with time-independent amplitudes 0„ and a master func- 
tion G{t) describing uniform oscillations of all the sites. 
After substitution of the ansatz (@J) into the Eqs. Q the 
following equation for the function G{t) is obtained: 

G + G = -CG^, (5) 
while the amplitudes (j)n satisfy algebraic equations: 

C0„ = -KY.^ [{C^n+l - ^nf + i^n-l - <f>nf] (6) 

where C is an arbitrary separation constant. Its absolute 
value can be always chosen to be equal to unity jl^l . 

While the dynamics of all the sites is governed by a 
unique function G(t), which can be easily found by inte- 
grating Eq. 10 , the spatial profile of possible solutions of 
Eqs. is determined by Eqs. © being of main interest 
for us. 



II. THE MODEL 

We consider a simple one-dimensional model of (non- 
linearly) coupled oscillators with the following Hamilto- 
nian: 

H = ^/i„^^|!| + yK) (1) 

where Un{t) is the displacement of nth unit mass oscil- 
lator from its equilibrium position, the constant s > 
characterizes the rate of spatial decay of long-range in- 
teractions between oscillators, and V{un) is given by 

? 4 

FK) = Y - f • 

The equation of motion for the displacement of the nth 
oscillator from its equilibrium reads: 
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We note that while the interaction decays algebraically 
for any finite power s, in the limit s oo we recover the 
case of short-range nearest neighbor interaction. 

The specifically chosen nonlinear potentials allow one 
to use the time-space separation technique 0, 0| , so 
that time-periodic solutions of (PJ can be written in the 
form 



III. BASIC TYPES OF COMPACT-LIKE 
DISCRETE BREATHERS 

It is important to note, that the dynamics of a DB, 
as well as its spatial profile, depend on the sign of G in 
Eqs. (|5lt)|l . This sign is fixed only for the case of un- 
coupled oscillators {K = 0) and for small values of K 
(C < 0), while generally it can be arbitrary. As a con- 
sequence, for large enough values of K Eqs. Q support 
two classes of DBs with different dynamical properties, 
since the sign of G defines the type of nonlinearity ('soft' 
or 'hard') in the Eq. ^ for the master function G{t). DB 
solutions of different classes possess essentially different 
core structures. As a consequence, the impact of long- 
range interactions on DB tail structure is also different 
for these two classes of DBs. 



A. Core structure of DBs 

In order to understand how two different classes of so- 
lutions of Eqs. Q appear, it is instructive to start with 
a simple case of three coupled oscillators. This simple 
model gives a rather good approximation for the sites of 
a DB core, which are practically not affected by the pres- 
ence of long-range interactions. Here we restrict ourselves 
in considering only symmetric DBs centered at site n = 
(it is straightforward to modify this approximate model 
for more complicated types of multi-site DBs). There- 
fore, by putting 0_i = 0i and (pi = kc^q, we finally ob- 
tain from Eqs. ((BJ the following expression for the central 
site amplitude: 



Un{t) = 0„G(t), 



(4) 



(7) 
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FIG. 1: Real roots of Eq. 0. The inset shows the correspond- 
ing values of 4>q/C computed from Eq. 0, squares, triangles 
and circles indicate the same quantity computed for the sys- 
tem of 201 coupled oscillators. 

while the coefficient k is a root of the fourth-order poly- 
nomial equation 

2Kk^ - {AK -f 1) 2,Kk^ + {K + 1)k - K = 0. (8) 

For each given real value of k obtained from Eq. the 
sign of C is fixed and determined by the condition of 
non-negative valued r.h.s. of Eq. iQ. 

In Fig. n] the real roots of Eq. JSJ are plotted in the 
range of the coupling constant K G [0,10]. There ex- 
ist two real non-trivial roots, whose absolute values stay 
below unity [and thus corresponding to single-site DBs 
with amplitudes of the central site (po greater than those 
of the neighboring sites 0±i] in the whole interval of non- 
negative values of K up to K ^ +oo. One of these roots 
originates from k = in the uncoupled limit (i.e. from a 
single-site excitation) and remains to be positive in the 
whole interval of K, see dashed line in Fig.^ The corre- 
sponding family of DBs have a non-staggered pattern of 
amplitudes - in-phase oscillations, see Fig. |21[a). The 
value of 0o/C computed for this root from Eq. ([TJl is neg- 
ative for < 2 (see dashed line in the inset in Fig. ^ , 
therefore one should choose C = — 1 for this type of DBs 
in Eqs. H5I6|I . Thus, DB solutions of Eq. (O with a non- 
staggered pattern should possess soft nonlinear proper- 
ties, i.e. their amplitude decays with growing frequency, 
as it follows from the Eq. ||SJ) with negative C. Here we 
mention, that the quantity (ji^/C computed for the 3-site 
model changes its sign above the value K = 2. However, 
the comparison with numerically obtained solutions of 
Eq. 10 for a larger system size {N = 201) indicates es- 
sential discrepancies when K ^ 1 (see triangles in the 
inset in Fig. . It comes from the fact that a DB core of 
breathers with the non-staggered profile extends while in- 
creasing K, approaching continuum compacton solutions 
[20j as ^ +00. Indeed, as shown in Fig. |21[a), while a 
DB core involves more and more sites with increasing 
its characteristic width in terms of the continuum coordi- 
nate X = n/K'^''^^ remains to be fixed and demonstrates 



rather good agreement with the value Lq « 2.92 reported 
for DBs in the continuum limit [20l| . 

Another non-trivial real-valued root of Eq. ijHJl origi- 
nates from K = — 1 in the uncoupled limit (i.e. from 
a staggered homogeneous excitation) and remains to be 
negative in the whole interval of K, see solid black line in 
Fig. n The corresponding DB family is characterized by 
a staggered profile in the core (while the tails have more 
complicated profile, as will be shown below), i.e. the cen- 
tral {n = 0) and the neighboring {n = ±1) sites oscillate 
in anti-phase. It is remarkable, that the corresponding 
value of 00 /C changes its sign at K = Kcr « 0.1, see 
solid black lines in the inset in Fig. ^ Therefore, for 
small enough values of the coupling constant K < Kcr 
both staggered and non-staggered types of DB solutions 
of (PI possess soft nonlinear properties, in accordance to 
the chosen type of the on-site nonlinear potential ||2Jl. 
The profile of the corresponding staggered-core DB so- 
lution at K ^ 0.07 is shown in Fig. I^^b). Notably, all 
the tail sites perform inphase oscillations in this type of 
DBs, similar to the case of nonstaggered DBs described 
above. However, the central site oscillates in anti-phase 
with all the rest of the lattice. 

The competition between on- and inter-site nonlin- 
earities results in the change of dynamical behavior of 
the DBs with staggered core profile as the coupling be- 
comes strong enough: The nonlinear term in Eq. Q 
for the master function G{t) becomes of the hard type 
for staggered DBs when K > Kcr- As the result, the 
amplitude of such a DB solution increases with grow- 
ing frequency. The corresponding profile is shown in 
Fig. |21[c). In this type of staggered DBs all sites in the 
core perform anti-phase oscillations. In the case of pure 
short-range interactions in the system, the staggered pat- 
tern persists for the whole spatial profile, including the 
breather tails. However, long-range interactions destroy 
the uniform staggered pattern, introducing a complicated 
domain-like structure in the breather tails, as will be ex- 
plained below. 

Unlike DBs with non-staggered profile, the staggered 
core DBs stay localized on a few sites as K increases, 
therefore the 3-site model gives a rather good approxi- 
mation for larger size systems for arbitrarily large values 
of K (see squares in the inset in Fig. Simultaneously, 
the amplitude of these DBs decreases as K ^ +oo (for 
a fixed value of the frequency). However, we note that 
there is no limit in frequency (and therefore, in energy) 
for this type of DBs, since the master function G{t) sat- 
isfies the Duffing equation (jSJ with the hard nonlinear 
term (C — 1). Note also, that as K ^ +oo, the in- 
fluence of the inter-site nonlinear interactions becomes 
more important than the effect of on-site nonlinearities 
for the oscillating in anti-phase sites of the DB core. 
Therefore, at large values of K the staggered-core DB 
asymptotically approaches the high-energy limit of dis- 
crete breather solutions in models with purely inter-site 
nonlinearities (Fermi-Pasta-Ulam lattices) 23]. 

Thus, for K > Kcr Eqs. (j2Jl support two different 
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FIG. 2: Profiles of various types of single-site (black color) and two-site (gray color) DB solutions of Eq. 0: (a) nonstaggered 
DBs corresponding to the branch marked with triangles in the inset in Fig. (S-type DBs), parameter values are: C = — 1, 
K — 1 (circles, squares) and K — 100 (triangles); (b) a staggered-core DB corresponding to the branch marked with circles in 
the inset in Fig. parameter values are: C — —1, K = 0.07 {K < Kcr); (c) staggered-core DBs corresponding to the branch 
marked with squares in the inset in Fig. (H-type DBs), parameter values are; C = 1, K = 1 [K > Kcr)- The long-range 
decay rate is s = 100, lines are provided to guide the eye. 



classes of DB solutions with soft and hard nonlinear 
properties, in what follows we will refer to such breather 
solutions as S-type and H-type DBs, respectively. The 
above discussed single-site DBs represent only particu- 
lar (basic) members of these two classes of solutions. In 
general, one can constract more complicated localized S- 
and H-type solutions - multi-site DBs. As an example, 
we mention here two-site DBs with two sites in the core 
oscillating with the same (maximum) amplitude, see gray 
lines and symbols in Fig. |21 The center of energy den- 
sity distribution is located in-between two sites in these 
DBs, so that they can be viewed as translated half site 
single-site DBs. 

Finally, we would like to remark, that the existence of 
the two different classes of DB solutions of Eqs. (|3Jl is 
the result of a competition between soft nonlinearity of 
the on-site potential and hard nonlinear inter-site inter- 
actions. Upon change the type of nonlinearity in the on- 
site potential [by changing the sign of the quartic term in 
Eq. ((SJ] , only the H-type DBs with staggered core profile 
survive. 



B. Linear stability 

Different dynamical properties of DBs with staggered 
and non-staggered core profiles result, in particular, in 
different stability properties of these excitations. In this 
section we perform linear stability analysis of basic types 
of H— and S— type DBs by studying the dynamical be- 
havior of a small perturbation e„(t) to a given DB solu- 
tion Unit). In order to construct a certain type of DB 
solution of Eq. Q, we solve numerically Eqs. © for the 
DB profile (t>n{t) with a certain sign of the separation 
constant C (chosen in accordance to the above performed 
analysis of DBs structure) and multiply it by a periodic 
solution G{t) of the Duffing equation according to 
the ansatz Q . Because of the time-space separation (0J , 
the stability properties remain qualitatively the same for 



the whole family of a given DB type with different fre- 
quencies f^B = 2tt/Tb [i.e. with different time-periodic 
functions G{t -t- Tb) — G{t)]. However, they might dras- 
tically change by varying the relative strength of the on- 
site and inter-site nonlinearities controlled by the cou- 
pling constant K. 

Thus, once a given DB solution u„(t) is obtained nu- 
merically, we add a small perturbation to it u„(<) = 
Unit) -\- en(t) and linearize equations of motion Eqs. (jS)) 
with respect to e„(t): 

e„ = ^ |(-u„+i - u„)^ (e„+; - e„) + (9) 

(it„_i - u„)^ (e„_; - e„)| - e„ + 3u^e„ . 

In order to simplify the stability analysis, here we re- 
strict ourselves in considering pure short-range interac- 
tion terms, thus we keep only the term with I = 1 in the 
r.h.s. of Eqs. which corresponds to the limit s — > cx). 
Being essential for DB tail characteristics, long-range in- 
teractions practically do not affect the core structure of a 
DB, provided that their decay rate s is sufficiently large 
(s > 3). Therefore, the impact of long-range interactions 
on stability properties of DBs is expected to be neghgible. 

The discrete breather acts as a parametic driver for 
small perturbations e„ (t) with the period T — Tb /2 be- 
ing the half-period of the DB solution [i.e. the half period 
of a given solution G{t) in Eq. Q]. 

Equations © define a map 

which maps the phase space of perturbations onto itself 
by integrating each point over the period T. Here we used 
the abbreviation x = {xi,X2, ...). The map H10|l is 

characterized by a symplectic Floquet matrix whose 
complex eigenvalues A and eigenvectors y provide infor- 
mation about the stability of the DB 2]. Here we note 
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FIG. 3: Evolution of the arguments of Floquet eigenvalues A 
with change of the coupling constant K for different types of 
single-site DBs: (a) S-type DB with the frequency fis = 0.9. 
Inset shows the Floquet spectrum ai K = 1, the unit circle 
in the complex plane is indicated to guide the eye; (b) H-type 
DB with the frequency Qb ~ 1.1. Inset shows evolution of 
absolute values |A|^ with change of K. A pair of eigenvalues 
corresponding to the unstable depinning mode (see the main 
body text for the details) of the two-site DB is indicated with 
gray color. 



that if all eigenvalues A are by modulus one, then the DB 
is linearly (marginally) stable. Otherwise perturbations 
exist which will grow in time (typically exponentially) 
and correspond to a linearly unstable DB. Upon chang- 
ing a control parameter (e.g. the coupling constant K) 
stable DBs can become unstable (and vice versa). Such 
a change of stability is appearing because two (or more) 
Floquet eigenvalues collide on the unit circle and depart 
from it 

In Fig. Ota) the typical Floquet spectrum is shown for 
a S-type single-site DB. All the eigenvalues A can be di- 
vided into three sub-categories: One pair of eigenvalues 
is always situated at exp(j7r) [denoted by squares in the 
inset in Fig. Ol^a)] ; it corresponds to perturbations along 
the DB periodic orbit (phase mode) and along the corre- 
sponding family of DB solutions 2] . The period of these 
perturbations coincides with the DB period. In addition, 
there are quasi-degenerated 28] bands of eigenvalues at 



A = exp(±iTB/2) [filled circles in the inset in Fig. (SJa)], 
which correspond to perturbations in breather tails. Such 
perturbations have characteristic frequency loi — 1 de- 
fined by the choice of the linear constant in the on-site 
potential l(5Jl, they would correspond to linear phonons 
if linear coupling between sites were introduced. Finally, 
there is a finite number of eigenvalue pairs bifurcating 
from the quasi-degenerated bands, which correspond to 
perturbations of the DB core sites [crosses in the inset 
in Fig. ISfa)]. The number of such isolated pairs is pro- 
portional to the characteristic DB core size, it grows as 
the coupling constant K increases, see Fig. ISJa). While 
increasing the coupling constant K, these isolated pairs 
move on the unit circle and new pairs bifurcate from 
the quasi-degenerated band, but they do not collide with 
each other and the non-staggered type of DB remains 
linearly stable up to the continuum limit K +oo. 

Notably, the Floquet spectrum of the S-type two-site 
DB is qualitatively the same as the one of the single-site 
DB for any K. The only principal difference is that it has 
two degenerated pairs of eigenvalues at A = exp(i7r), since 
in the uncoupled limit A' = the corresponding solution 
has two sites excited with equal amplitude. Usually such 
degeneracy of eigenvalues is lifted for any nonzero K , and 
one pair of eigenvalues is "pushed out" from A — exp(z7r) 
either along the real axis or along the unit circle. How- 
ever, in the case of purely nonlinear interactions between 
sites the symmetric two-site DB "cuts" the effective lin- 
ear chain ^ into two non-interacting halves. Therefore, 
the additional degeneracy of eigenvalues corresponding 
to symmetric and anti-symmetric perturbations with re- 
spect to the DB center remains for any K. 

In contrast to the S-type DBs, the H-type single-site 
and two-site DBs are linearly stable only within cer- 
tain windows of the coupling constant K values, see 
Fig. |2Ib) . Close to the critical value of the couphng con- 
stant K = Kcr, below which the H-type DBs do not exist, 
both single- and two-site H-type DBs experience strong 
instabilities connected with tangent bifurcations of these 
solutions with other ones, having more complicated spa- 
tial structure. In addition, there is another instability of 
a finite strength, appearing in certain windows of the pa- 
rameter K, see inset in Fig. I^Jb). In general, apart from 
several small intervals in K, for any given value of K 
only one of these two DB configurations is stable. The 
corresponding unstable perturbation - the "depinning" 
mode - "tilts" the single-site (two-site) DB towards the 
half-site shifted stable two-site (single-site) one. Chang- 
ing the coupling constant, the stable configuration varies 
from the two-site to single-site DB and back [at iiT « 0.15 
and K w 0.3, see Fig. I^Jb)], so that the exchange of sta- 
bility process '25^ is observed. This exchange of stability 
process can be connected to an exchange of the domi- 
nant roles between inter-site and on-site nonlinearities. 
Indeed, typically for models with purely inter-site non- 
linearities (Fermi-Pasta-Ulam latti ces) the basic stable 
configuration is the two-site DB Q,|2^. In contrast, for 
models with weak coupling between sites and a non- linear 
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onsite potential in the form j^J (Klein-Gordon lattices) 
the single-site DB configuration is stable 2]. 

Of principal interest is an influence of the unstable 
depinning mode on the dynamical behavior of staggered- 
core DBs. A small perturbation along this mode is gener- 
ally known to result in depinning of the unstable DB from 
its initial position. Depending on the relative Hamilto- 
nian energy ^ of the perturbation (AH/Hdb) and on 
the strength of instability of the depinning mode, the re- 
sulting behavior might vary from quasiperiodic-like oscil- 
lations between two neighboring unstable positions (in a 
well of the corresponding " Peierls-Nabarro potential" ) to 
quasi-regular or even chaotic-like motion along the chain 
(see e.g. |2^|2^). Generally, the depinned DB resonates 
with linear phonons through its excited internal modes 
and starts to radiate energy. Therefore, eventually it will 
be trapped again at some stable position or even disap- 
pear completely transferring totally its energy to excited 
delocalized phonons. 

However, in the case of purely nonlinear dispersion 
there are no linear phonons in the system, and all pos- 
sible linear resonances are suppressed. As a result, one 
can observe almost perfect quasi-periodic oscillations of 
perturbed unstable DBs between two neighboring stable 
positions, see Fig. ^a). Similar quasi-periodic behavior 
is observed in dynamics of a stable DB with perturbation 
along one of its internal modes, see Fig.0Jb,c). In the lat- 
ter case no detectable radiation (within the used double 
precision) of energy from the DB core was observed dur- 
ing dynamical simulation over 10^ breather periods. In 
this respect, a one-dimensional Hamiltonian lattice with 
purely anharmonic interactions between sites ^ might 
be an interesting " toy model" to study more complicated 
exact DB solutions like quasi-periodic and moving dis- 
crete breathers. 



IV. LONG-RANGE INTERACTIONS EFFECT 

Let us now fix the value of A' = 1 and study the in- 
fiuence of long-range interactions on the spatial profile 
of a DB. In Fig. [S] the profiles of the H-type single-site 
DBs are shown for various values of the decay constant s. 
They were obtained by solving Eqs. O numericallv with 
use of the standard Newton scheme |2J| for the chain of 
A'' = 201 oscillators (-100 < n < 100) with periodic 
boundary conditions [23. Similarly to the case of non- 
linear short-range interaction, the computed DBs have 
a compact-like structure with mainly three central sites 
oscillating (see Fig. (SJ. The other oscillators are almost 
at rest and can be considered as breather tails. The pres- 
ence of long-range interactions breaks the uniform super- 
exponential law of the spatial tail decay known for the 
case of pure short-range nonlinear interactions [3 ' 
introducing several cross-over lengths. 

For a few central sites, lying within the breather core, 
the forces due to the long-range interactions are negligible 
as compared to those due to the nearest-neighbor interac- 



tions. Thus, the central part of a breather is practically 
not affected by the presence of long-range interactions. 
However, at some distance Li from the DB center inter- 
actions with the nearest neighbors (having small enough 
amplitudes) become of the same order as the long-range 
interactions with the DB core (central three sites having 
the highest amplitudes). This distance is the first cross- 
over length, where the long-range interactions come into 
play. It can be roughly estimated by an assumption, 
that at the distance Li interactions with the breather 
core sites are exactly compensated by interactions with 
the nearest neighbors. Thus, keeping only the leading 
order terms in the sum in the r.h.s. of © one obtains: 



(11) 



Since for \n\ < Li the relations between the ampli- 
tudes (j)n are practically the same as in the case of pure 
nearest-neighbor interactions, i.e. they follow the super- 
exponential law 0„ ~ (_l)"|<^^|exp(ln(3)|n-l|)^;^ < |„| < 

Li, one can obtain from Hll|l : 



|0^|cxp[31n3(Li 



-2)l-3(Li - 1)^ 



1. 



(12) 



In the limit of extremely large values of s the distance Li 
will be also large, and satisfy 



Li 



In s 

hTa 



(13) 



Thus the first cross-over length Li grows approximately 
logarithmically with s. The numerical results in Fig. [S] 
yield Li{s = 10) w 3, Li(s = 20) « 4, Li(s = 50) « 5 
and Li{s = 100) ~ 6. They compare very well with the 
corresponding solutions of lO: 2.71, 3.46, 4.44, 5.17. 
Therefore, even extremely fast (but still algebraically) 
decaying in space long-range interactions essentially de- 
stroy the concept of compact-like breathers, since only 
amplitudes of a few sites in the tails obey the super- 
exponential law of decay, while the rest of the tail am- 
plitudes decay much slower in space (following a power 
law, as will be shown below). 

At large distances from the breather center n ^ 1 (due 
to the single-site DB symmetry around n — we consider 
here only non-negative values of n) the impact of short- 
and long-range interactions is exchanged: now the most 
powerful contribution comes from the interaction with 
the breather core, while nearest neighbors, due to their 
small amplitudes, practically do not affect the dynamics 
of a tail site. Thus, for large n one can derive from ® 
the following asymptote: 



K 

"c 



{n-iy (n-l-1)* 



(14) 



which in fact gives a rather good approximation for all 
tail sites starting from the first cross-over point n = Li 
(see solid lines in Fig. E)). Note that only amplitudes of 
the two DB core sites and the sign of separation constant 
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(c) 



FIG. 4: Dynamics of a perturbed DB: (a) Position of the DB center Xc, calculated through the energy dencity /i„ 
Xc ~ '^^{n ■ hn) / "^^hn. At time t = G the single-site H-type DB {Q,b = 1.1, K — 0.3) was perturbed along the unstable 
depinning mode with the relative energy of perturbation AH/Hdb ~ 10"*; (b) Energy density h„ Q evolution in dynamics 
of the stable single-site H-type DB {Qb = 1.1, K = 0.2) with a symmetric perturbation along one of the DB internal modes. 
The relative energy of perturbation AH/Hdb ~ 10~^; (c) The Fourier transform of Mi2(t) (the DB is centered at n = 11) for 
the same dynamical simulation as in (b). Vertical dashed lines indicates the first and the third DB harmonics ^Ib = 1.1 and 
3Qb = 3.3. 
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5 10 15 20 25 30 



n 

FIG. 5; Single-site H-type DB solutions of with C = 1 and 
K — 1 for various exponents of the long-range interaction: 
s = 10 (squares), s = 20 (stars), s = 50 (crosses) and s — 100 
(triangles). 



C are needed to obtain this asymptote for tail amplitude 
distribution. In this respect we found the simple three- 
site model, discussed in Sec. Illll to be very fruitful: it 
gives full information not only about the DB core sites, 
but about tail characteristics as well. 

Note that the specific structure of a staggered DB 
core with a central site n = and two neighboring sites 
n = ±1 having amplitudes 4>n of opposite signs stipulates 
several other cross-over lengths connected to changes of 
the sign of the r.h.s. in H14|l which manifest as sin- 
gularities in the logarithmic plots in Fig|^ The most 
pronounced cross-over a,t n = L2 is associated with the 
change from a single power law to a more complex 
one ifTH) . see Fig. El Indeed, in the case n 3> 1, (s/n) <C 1 
the expression p4fl can be re-written as: 

*..-M{l + 2„= + 4 + o(i;)}, (15) 




FIG. 6: Single-site H-type DB solutions of with C — 1 
for various exponents of the long-range interaction: (a) s — 
100; (b) s = 20. Circles: numerical results. Solid lines: tail 
asymptotes 1141 . Dashed lines: location of Li. The inset in 
(b) indicates the change of amplitudes sign around the cross- 
over point L2. 



(C = 1 , K < for a staggered core DB and C = — 1 , k > 
for a non-staggered DB). Thus, in leading order, at large 
enough distance from the DB center n > L2 the tail 
amplitudes follow the same power law n"" as the decay 
of long-range interactions. Since L2 is defined by the 
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vanishing of the bracket on the right hand side of Eq. 1151) 
we obtain in leading order 

The corresponding values of L2 for s — 10, 20, 50, 100 and 
K = -0.382 with ^ are 5.6,11.2,28,56. They com- 
pare reasonably well with the numerically observed ones 
L2(s = 10) « 7, L2{s = 20) « 13, L2{s = 50) « 31 
and L2{s = 100) « 64. In-between the two characteris- 
tic length scales Li < n < L2 the tail amplitudes decay 
following a more complicated power law p4|l. 

V. CONCLUSIONS 

To conclude, we revealed the influence of long-range 
nonlinear interactions on the spatial profile and proper- 
ties of compact-like discrete breathers in a model of cou- 
pled oscillators with pure nonlinear dispersion. As we 
demonstrate, it is the intriguing property of the model 
under consideration, that it supports two classes of dis- 
crete breathers - with staggered and non-staggered spatial 
profiles of a DB core - having different dynamical proper- 
ties. The dynamics of a non-staggered DB is essentially 
governed by the soft nonlinear on-site potential, while 
DBs with the staggered core have the opposite, hard, 
type of nonlinear dynamical behavior caused mainly by 
the presence of nonlinear interactions in the chain. Apart 
from different dynamical properties, the influence of long- 
range interactions on spatial profiles of these two types 
of DBs is also different. With the algebraic spatial de- 
cay long-range interactions introduce a new length scale 
which becomes essential at large enough distances from 
a DB core. We show, that the effect of long- and short- 
range terms competition results in the appearance of a 
characteristic cross-over length Li in both types of DBs, 
at which the spatial tail decay drastically changes from 
the super-exponential law to the algebraic one. For large 



powers s of the long-range interactions spatial decay the 
cross-over length Li scales logarithmically with s. The 
tail asymptote demonstrates complex power law spa- 
tial decay, which follows essentially the same algebraic 
decay as the long-range interactions in the system at 
large enough distances from the DB core. While for 
non-staggered DBs the influence of long-range interac- 
tions manifests through the only characteristic cross-over 
length Li, the spatial profile of DBs with the staggered 
core possess several other cross-over lengths associated 
with sign changes of the asymptote H14|l . Thus, the spa- 
tial pattern of oscillations in a DB with the staggered 
core becomes rather complicated in the presence of long- 
range interactions: its core sites perform anti-phase os- 
cillations, while its tails are splitted in several domains 
of in-phase oscillations. 

Finally, we would like to mention that the discussed 
case of purely nonlinear coupled oscillators represents a 
simple model to reveal properties of nonlinear excitations 
in a system without linear phonons. As we demonstrated 
in this paper, "switching off' the phonons leads not only 
to the change in characteristic rate of spatial localization 
of energy, but to appearance of several other intriguing 
dynamical properties of nonlinear localized excitations. 
Especially it eliminates the possible source of linear res- 
onances which otherwise would destroy quasiperiodic 
breathers and possibly also moving breathers. The dis- 
cussed models thus allow to obtain a better understand- 
ing of the general problem of existence/non-existence of 
quasi-periodic and moving discrete breather solutions. 
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